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The highly spotted photosphere of the young rapid rotator Speedy 
Mic 
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ABSTRACT 

We present high resolution images of the young rapidly rotating K3 dwarf Speedy Mic 
(BO Mic, HD 197890). The photospheric spot maps reveal a heavily and uniformly spotted 
surface from equatorial to high latitude regions. Contrary to many images of similar objects, 
Speedy Mic does not possess a uniform filling at high latitudes, but exhibits structure in the 
polar regions showing greatest concentration in a particular longitude range. 

The asymmetric rotation profile of Speedy Mic indicates the presence of a companion 
or nearby star which shows radial velocity shifts over a timescale of several years. Using a 
simple dynamical argument, we show that Speedy Mic is unlikely to be a binary system, and 
conclude that the feature must be the result of a chance alignment with a background binary. 

Complete phase coverage on two consecutive nights in addition to 60% phase coverage 
after a three night gap has enabled us to track the evolution of spots with time. By incorporat- 
ing a solar-like differential rotation model into the image reconstruction process, we find that 
the equator laps the polar regions once every 191 ± 17 d. This finding is in close agreement 
with measurements for other late-type rapid rotators. 

Key words: Line: profiles - Methods: data analysis - Techniques: miscellaneous - Star: 
Speedy Mic (HD 197890) - Stars: activity - Stars: atmospheres - Stars: late-type 



1 INTRODUCTION 

The distribution of starspots on young rapidly rotati ng star s has 
been determined for a number of objects (Strassmeier 200*2); and 
in some instances we now have photospheric image reconstructions 
at several epochs on a single object. While patterns in the distribu- 
tion of emergent flux, as a function of spectral type for example, are 
just becoming discernible, it is still not entirely clear how the dis- 
tribution may change on a single object with time. Direct evidence 
for chromospheric variation on timescales similar to the solar mag- 
netic activity cycle have shown that some moderatel y rotating stars 
possess similar cycles of th eir own <Wilsorj[l978l iBaliunas et"ail 
1 19951 iDonahue eydj^9^). Similarly, period variations derived 
from photometry <Hallll99lblHenrv et aljl995lllvlessina & Guinart 
2002) have indicated that rapid rotators also show periodic varia- 
tion in the emergent latitude of starspots. While long term variation 
in photometric brightness indicates cyclic activity, no evidence for 
variation of emergent latitudes in Doppler images has yet been pre- 
sented, since very few objects have been observed consistently over 
a long enough time basis. 

Here we present new observations of one of the fastest and 
bright est young ra pidly rotating stars in the solar neighbourhood. 
Mont e~t alJ 1200 ll) have shown that Speedy Mic, like its rapidly 
rotating cousins AB Dor and PZ Tel is indeed a member of the 



E-mail: jrb3@st-andrews.ac.uk 



Local Association. Speedy Mic (HD197890) is a K3 dwarf star 
with an equatorial rotation velocity, v sin i ~ 128 km s" 1 . We first 
observed this star in 1998 jBarnes et alll200ll) (B01), finding that 
spots are distributed at a range of latitudes. In B01 we gave a sum- 
mary of various observations of Speedy Mic. Since the first de- 
tectio n of magnetic activity iBromag e et alj|l992llMatthews et all 
1 1994) . it has been shown that S peedy Mic is one of the most 
active solar neighbourhood stars ISingh et al]|l999 |). with a n X- 
ray luminosity of Log(L x / Lbol) = -3.07. ISingh et alJ 1 19991) also 
found evidence for flaring activity and significantly subsolar coro- 
nal Fe abundances in accordance with obs ervations of oth er ac- 
tive stars. A more recent measurement by (Mak arovlEool indi- 
cated an X-ray luminosity of L x = 116.78 x 10 29 ergs s" 1 , an 
order of magnitude higher than that of the well studied K0 star 
AB Dor (L x = 15.22 x 10 29 ergs s" 1 ) for example. This implies 
hog(L x I Lbol) ~ -2, which is an order of magnitude higher than 
measurements for similar stars, again indicating probable flaring 
activity at the time of observations. 

In B01 we also confir med the presence of Hq transients, first 
reported bv ljeffrie"sl 1 19931) . These transients are thought to be the 
result of clouds of cool material, analogous to solar prominences, 
passing at large distances from the rotation axis in front of the stel- 
lar disc. Both studies found prominences at or above the co-rotation 
radius of the star. The Ha profile in the present data set is studied 
in a separate publication ( Dunstone et al. 2005), revealing a densely 
packed prominence system beyond the co-rotation radius. 
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Table 1. Observation log for 2002 July 18, 19 & 23 



Object 


UT Start 


UT End 


Exp time 
[s] 


No. of frames 


comments 


2001 Jul 18 


HD106911 
HDl 26053 
G1673 
G1472 
G1701 
Speedy Mic 
G1849 
HD 16160 


08:41 
08:46 
08:53 
09:01 
09:57 
10:10 
19:45 
20:01 


08:42 
08:51 
08:58 
09:06 
10:07 
19:42 
19:58 
10:06 


40 
250 
300 
300 
600 
300 
800 
300 


1 
1 
1 

102 


B5V telluric standard 
G1V RV standard 
K5V spectral standard 
K1V spectral standard 
M1V spectral standard 
Target 

M3.5V spectral standard 
K3V spectral standard 


2002 Jul 19 


HD106911 
HDl 19850 
Speedy Mic 
G1825 
HD4628 


08:07 
08:14 
10:05 
19:19 
19:26 


08:08 
08:24 
19:18 
19:24 
19:29 


60 
600 
300 
300 
180 


1 
1 

97 
1 
1 


B5V telluric standard 
M1.5V spectral standard 
Target 

K7V spectral standard 
K2V spectral standard 


2002 Jul 23 


HD106911 
Speedy Mic 
G1825 
HD4628 


08:02 
09:46 
19:19 
19:26 


08:03 
18:54 
19:24 
19:29 


40 
300 
80 
180 


1 

97 
1 
1 


B5V telluric standard 
Target 

K7V spectral standard 
K2V spectral standard 



In this paper we present images at closely separated epochs. 
We examine the reliability of the strength of reconstructed features 
in the maps through use of a bootstrapping method. This allows 
us to distinguish between real evolution of starspot features and 
uncertainties in the reconstru ction process . We c ompare our im- 
ages with those published by IWolter et all J20051) (WSW05), de- 
rived from observations made between two and three weeks after 
those presented here. We also combine the complete timeseries and 
determine the differential rotation rate at the photospheric level. 



2 OBSERVATIONS 

Observations were made with the University College London 
Echelle Spectrograph (UCLES) at the Anglo Australian Telescope 
(AAT) on 2002 July 18, 19 & 23. With a 1.2" slit width, a mean 
spectral resolution of ~ 44000 at 5450A was attained with cover- 
age from 4376 A to 6892 A in 47 orders on the 2Kx4K EEV CCD. 
This corresponds to ~ 6.7 km s _1 at the mean wavelength. The 
pixels were binned by a factor of two in the wavelength direction 
giving two pixels per spectral resolution element. Exposure times 
for Speedy Mic were 300 sees for all spectra. The seeing varied 
between 0.9" and 2" with a mean of 1.5". Cloud gradually built 
up to 80% cover during the second half of July 23. This resulted 
in loss of target counts and contamination of spectra by moonlight. 
We have corrected for t his in a similar manner to that des cribed in 
iDonati et allfcOOOl) and barnes & Collier Cameror]<200ll) . We do 
not detect any further significant changes in continuum level during 
the course of observations, which may be due to flaring activity for 
instance. 

We observed a number of slowly rotating template stars with 
effective temperatures similar to the photosphere and spots of 
Speedy Mic. These spectral standards are used to represent the local 



intensity profiles in the image reconstruction process. A featureless 
rapidly rotating B5V star (HD 106911) was observed as a telluric 
standard and also served the purpose of tracing the spectral orders 
in the extraction process. 

3 DATA REDUCTION 
3.1 Extraction 

Pixel to pixel variations were removed using flat-field exposures 
taken with an internal tungsten reference lamp. The worst cosmic 
ray events were removed at the pre-ex t raction stage using the 
FIGARO routine B CLEAN IShortridgel [l993). Scattered light 
was modelled by fitting polynomials of degree 7 to the sets of 
inter-order pixels at each X-position in each frame. The spectra 
were extracte d usin g ECHOMOP, the echelle reduction package 
developed by Mills ( 1994). The Thorium-Argon arc-frames used 
for wavelength calibration were extracted in conjunction with a 
target spectrum, and calibrated using this package. The orders 
were extracted using ECHOMOP's implemen tation of the optimal 
extraction algorithm developed by Hornj <1986h . ECHOMOP 
propagates error information based on photon statistics and readout 
noise throughout the extraction process. 



4 THE SPECTRA 

We applied leas t squares deconvolution jDonati et alJ Il997t 
iBarnes et al J 19981) to each spectrum, thereby deriving a single pro- 
file at each observation phase. The list of lines used for deconvolu- 
tion w as taken from the Vienna Atomic Line Database l Kup ka et"al] 
1 19991) for a star with T = 4750 K and log g = 4.5. In total, 5680 lines 
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Figure 1. Time series spectra for 2002 July 18, 19 & 23, with phase plotted against velocity. The mean profile has been subtracted from each time series. 
White features correspond to starspot signatures. 



Table 2. Mean S/N statistics for each night of observations. Cloud was re- 
sponsible for poorer S/N ratio in the latter part of the night of 2002 July 23. 
We thus give statistics in two parts for this night. 



Table 3. System parameters for Speedy Mic 



Date 


No. of spectra 


Input S/N ratio 


Output S/N ratio 


Gain 


18 


100 


77.0 ±5.5 


3729 ± 493 


48 


19 


97 


74.3 ±11.5 


3692 ± 766 


50 


23 


51 


86.8 ± 4.4 


4211 ±438 


49 


23 


16 


43.7 ± 14.8 


1784 ± 969 


41 



with depths between 0.05 and 1.0 were used in the deconvolution. 
Due to overlap of spectral coverage of adjacent orders in the blue 
part of the spectrum, a total of 7547 images were used to derive 
the rotation profile. The noise statistics are presented in Table [2] 
The phased timeseries spectra are plotted in Fig.Qwhere the mean 
profile has been subtracted to reveal many starspot (white trails) 
signatures. 



4.1 System parameters 

The system parameters for Speedy Mic have been determined by 
a number of authors, including B01, in a previous paper based 
upon spectra taken during 1998 July. During the 1998 observations, 
only a few spectra at a number of different instrumental config- 
urations (at two telescopes) were available. However the data set 
presented here provides us with nearly three complete rotation cy- 
cles with very high quality data. Using our Doppler imaging code, 
we combined all 264 usable spectra into a single timeseries and 
re-determined the radial velocity, v ra d, inclination, i, equatorial ro- 



P[d] 

v r [kms -1 ] 
v sin i [kms -1 ] 
Axial inclination [deg] 



0.38007 ± 0.00005 
-8.0 ± 1.0 
132 ±2 
70.0 ± 5 



tation velocity, v sin i, and period, P. The period measurement does 
not include differential rotation shear which is measured and dis- 
cussed in a later section. We also emphasise that the period as mea- 
sured relies on the repetition of the starspot trails in the timeseries, 
from night to night, and thus represents the mean period of Speedy 
Mic. The quoted uncertainty is for the measurement on this data 
set and does not take into account any changes which may occur if 
the mean rotation rate changes. In other words, if a magnetic cycle 
is present then the mean location of spots may change leading to 
a different rotation period measurement. The period is nevertheless 
in agreement with the 0. 380 ± 0.004 d photometric period found 
bv lCutispoto et alJil99% 

As with the 1998 data set, we deconvolved the spectrum of the 
slowly rotating K3 star HD 16160 in the same manner as our target 
spectra. This profile was used to represent the local intensity profile 
in the imaging process, and was scaled according to foreshortening 
cosine and limb darkening. Wc used a 4-paramctcr limb darkening 
law jciaret 2000) interpolated to a temperature of T = 4890 K. 

The mean profile from all spectra indicates a significant asym- 
metry in the rotation profile in that the red wing of the profile ex- 
hibits excess absorption. We discuss this feature at length in j5|in 
light of this finding and previous high resolution spectroscopic ob- 
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Figure 2. Mean profile from all spectra. The 1998 July 11-13 profile is 
plotted with the mean profile derived from observations in 2002 July 18, 19 
&23. 



servations made on Speedy Mic. Since we are certain that this ab- 
sorption feature does not arise from Speedy Mic itself, we have 
excluded it from the imaging process by artificially inflating the er- 
ror bars (effectively to infinite size) in the region of 85 km s" 1 to 
150 km s" 1 . The system parameters are presented in Table [3] and 
indicate values which are close to our original estimations. Notably, 
the preferred inclination is i = 70° rather than the i = 55° which we 
initially measured. The v sin i measurement is 3% greater than that 
measured on the 1998 data set. We attribute this to removal of the 
redshifted absorption feature which effectively removes the right 
hand portion of the profile. This includes the region of strong cur- 
varture where the line meets the continuum, normally providing a 
strong constraint on the extent and therefore v sini of Speedy Mic. 
Since we had to perform a similar task on the blue wing of the 1998 
data set, where the absorption feature appeared even more signifi- 
cant, it is likely that there is some subjective element in choosing 
which region to exclude from the fit. If too much is excluded, we 
are likely to slightly underestimate the v sin i while excluding too 
little may overestimate v sin i. Given the higher and more consistent 
quality of the current data set we will assume v sini = 132 km s - 
in this paper. The 3% difference in ?;sini will nevertheless have 
little effect on the image reconstructions. 



5 SPEEDY MIC: SINGLE STAR OR BINARY? 
5.1 Real or artifact? 

The mean profile in our timeseries indicates excess absorption in 
the red wing. A similar feature was present in the 1998 data (B01) 
which we attributed to the presence of a background star in the slit 
at the time of observations. The re-appearance of this feature in 
the red wing of the rotation profile (Fig. [2} suggests the possible 
presence of a companion star, indicating that Speedy Mic may in 
fact be a binary or multiple system. In order to ascertain whether 
the absorption feature is real and not an artifact, or a problem with 
the deconvolution process, we deconvolved various combinations 
of blue and red orders of each spectrum independently for the night 
of July 18. In all cases, the excess absorption in the red wing is 
present. We also deconvolved the full spectrum using a slowly ro- 
tating template star rather than a line list in case mis-matches in 
the synthe tic line l ist we re causing problems. This method, de- 
scribed in Bame^ 120041) . uses singular value decomposition to 



perform the deconvolution. Here, it is possible to obtain a much 
closer match between the reconstructed spectrum (i.e. template re- 
convolved with rotation profile) and the observed spectrum. Apply- 
ing this method of deconvolution nevertheless also resulted in the 
appearance of the absorption feature. 



5.2 Nearby objects 

Having confirmed that the absorption feature is real, we must con- 
sider that a nearby star was present in the slit of the spectrograph 
during the observations, or that Speedy Mic is in fact a binary or 
multiple system. We have searched the 2 Micron All Sky Survey 
(2MASS) 1 catalogue which does not list any objects close enough 
to Speedy Mic which could contaminate our spectrum. The 2MASS 
J-band image shows an object (unidentified by the survey) at a dis- 
tance of approximately 11" and position angle of 13°. Compared 
with objects in the same field, it appears to possess a magnitude 
of J ~ 13-14 and thus it seems unlikely that even if it is a binary 
system, it would contaminate our spectra. Similarly, the seemingly 
redder object which appears only in the K-band (K = 13 - 14) image 
at 14" and position angle 185° is unlikely to have appeared in the 
slit. 

A s urvey for brown dwarf co mpanions to nearby stars has been 
made by McCarthy & Zuckerman 1 2004) using coronagraphs to oc- 
cult the primary star. Below 5" the limit of detectability decreased 
almost linearly with distance, reaching a sensitivity of ~ AJ = 2 
at 0.15" separation. Beyond 5", the limiting AJ = 12.5 remained 
approximately constant, however no candidate brown dwarfs were 
detected close to Speedy Mic in the survey. 



5.3 Is Speedy Mic a binary? 

Observations of Speedy Mic in the past have attributed it single 
star status, nevertheless a range of rotation velocities have been re- 
ported. From the FWHM of the Ca II K line. lBromage et alJ l 1992) 
estimated v sini = 120 ± 40 km s -1 and v ra d = -6.5 km s _1 with 
no systematic variations in the latte r value over the six days dur- 
ing which observations were made. I Anders et alJ 1 19931) however 
reported v sini = 170 km s~ using spectral synthesis of two lines 
centred around 6710 A. By combining data from simultaneous 
data at two observatories, B01 found that v sini = 128 kms -1 and 
v ra d = -8.0 km s _1 . A range of v rac i of -8.0 to -10.5 for individ- 
ual nights of observations and instrumental setups was most likely 
due to incomplete phase coverage and inaccurate internal calibra- 
tions. The estimated error on this latter measurement is of order 
2 km s _1 . As was shown in £14,11 we found v sini = 132 km s~ 
and v ra d = -8.0 km s _1 . 

We are therefore presented with the contradictory observations 
of no (or at most small) radial velocity variations and an absorption 
feature, apparently due to a secondary component, which appears 
in the blue wing at one epoch, and the red wing at another. To assess 
whether Speedy Mic could plausibly have a low mass companion, 
we have estimated the expected radial velocity amplitude of each 
component for a range of secondary stellar masses. With an as- 
sumed binary orbital inclination of 70° and a circular obit, the mass 
function is defined <Hilditchl200ih as 



rri2siri 3 (i) 
(m x + m 2 ) 2 



1.0361 x IQ^KfP. 



(1) 



1 See http://irsa.ipac.caltech.edu 
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where masses are given in M®, period P in days and Ki 
in km s _1 . Table [4] gives K\ velocities for a range of secondary 
masses and periods. Since the absorption feature is not seen to 
move during the 5 night span of observations, we have assumed 
that the proposed star is at quadrature, and estimate that it can- 
not have moved through more than 10% of the orbit. This yields a 



lower limit to the orbital period of 50 days. The maximum period 
of 8 years is assumed from the approximation that the absorption 
feature took 4 years (1998 to 2002) to traverse from the blue wing 
(B01) to the red wing (present observations). 

Based on measurements made by various authors and our own 
measurements in (B01), we estimate an upper limit of 4 km s _1 to 
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Figure 4. Standard deviation maps from bootstrap reconstructions of 100 data sets. Note that the range of plotted values is to 0.25 (plotted for clarity), so the 
greyscale is not directly comparable with that used in Fig. [3] 



the radial velocity variations exhibited by Speedy Mic. This con- 
strains an upper spectral type for the secondary component of M5 in 
the considered period range. In order that the secondary component 
exhibits sufficient radial velocity variation to explain the observed 
absorption features, we can narrow down the secondary spectral 



type to late M and a binary orbit of tens of days (the M8, P=50 d 
case in Table|4j. The proposed secondary must itself be a rapid ro- 
tator and not tidally locked with the primary. However, the binarity 
argument fails because an M8/K5 pair would exhibit a large lumi- 
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Table 4. K± and K2 velocities for a range of secondary masses and orbital 
periods of 50 d and ~8 yrs. The primary mass is ass umed to be 0.8 Mq. 
Secondary masses and spectral types estimated from <Coxl2000l) . 



Secondary 




Kx 


K 2 


Sp Type 


Mass [M & ] 




P [days] 




in, 2 


50 


3000 


50 3000 


MO 


0.51 


21.1 


5.91 


36.2 9.3 


M2 


0.40 


19.2 


4.91 


38.5 9.8 


M5 


0.21 


11.3 


2.89 


43.1 11.0 


M8 


0.06 


3.6 


0.92 


48.0 12.2 



nosity ratio, with log(Ls pe edy/LMs) ~ —3.3 iGolimowski et alJ 
2004), ensuring that the secondary profile would not be seen. 

The only remaining possibility is that Speedy Mic is in a 
chance alignment with a binary. We find no evidence for starspot 
transients in this part of the profile, although they may be below 
the detection limit. We therefore conclude that the additional ab- 
sorption feature arises from a background object which has re- 
mained undetected in present surveys due to its close proximity 
to Speedy Mic. 



6 THE IMAGES 

The reconstructed images of Speedy Mic are presented in Fig.[3]for 
each individual night of observations. All images are phased on the 
epoch HJD = 2452473.5, midnight on 2005 July 18. As expected 
from the numerous trails in Fig. Q many starspot features are re- 
constructed at most latitudes and are found uniformly spread at all 
phases. The images from the nights of July 18 and 19 show many 
similarities, where the same spots are seen 2-3 rotations apart. By 
July 23 most spots at low and intermediate latitudes are also still 
visible whereas there may have been some evolution at higher lati- 
tudes. 

6.1 Image reliability 

We have attempted to obtain some measure of the reliability of the 
reconstructed images by applying a Monte-Carlo bootstrap proce- 
dure fEfron 1997t |Press et alJ 19921) to each night of observations. 
This involves generating a large number of data sets from our orig- 
inal data set and reconstructing images as follows. Firstly, we pre- 
serve the order of the rotation profiles. For the n profiles, in a data 
set, we pick at random and with replacement n profiles. If a given 
profile is chosen m times, we divide the errors on that profile in the 
new data set by \/m, whereas profiles not chosen are given very 
large error bars in order to exclude them from the timeseries. In this 
manner, we generated 100 data sets for each night of observations 
and reconstructed an image. This technique therefore requires no 
manipulation of the original profiles, and involves re-distributing 
the error weights only. As such, the amount of noise in the data is 
conserved. We can thus fit each data set to the same level of \ 
in order to obtain a reliable estimate of the uncertainties in recon- 
structed pixel values in the maps. 

An image with values representing the standard deviation 
for each pixel was generated and is shown in Fig. [4] It is clear 
from these maps that most of the strong features are reliably re- 
constructed as the cores of the features show uncertainty in filling 
factor of typically 0.1. As expected, the elongated features at low 



latitude which typically show a spot occupancy of 0.1 also show 
high uncertainties in spot filling. These features are a reflection of 
the uncertainties in the rate at which starspots pass through the pro- 
file. The gradient is proportional to the sine of the latitude, 9, and 
sin(8) is a rapidly changing function near the equator. 

The main point of this exercise is to distinguish between evo- 
lution of spots and uncertainties in the reconstruction process. Fo- 
cusing on the four spots seen at phase 0.8 to 0.9 and centred around 
latitude 60° in the July 18, we see that the arc-like structure joining 
the spots possess a greater uncertainty than the spots themselves. 
By July 19, there does appear to have been genuine evolution of 
these features, especially the two spots at phase at phase 0.8, al- 
though the reliability of the lower latitude spot of these two features 
is quite large on July 19. It is nevertheless clear that there are many 
features common to all images which makes us confident that we 
can reliably measure the differential rotation. 

6.2 Differential rotation 

The large number of repeated features and high effective sur- 
face resolution (full width at half maximum of mapping profile is 
7.69 km s -1 , corresponding to a resolution of 5.2° on the equa- 
tor) make Speedy Mic an ideal target to attempt a measurement 
of differential rotation. By optimising the degree of fit of a solar- 
like differential rotation law (where Q(9) = Qq — Affein (9), 
6 is the stellar latitude and fin the equatorial rotation ve- 
locity) to the full data set, incorporating the full time series 
of data, we find that fi = 16.5361 ±0.0006 rad.day" 1 and 
An = 0.03295 ± 0.00294 . This equates to an equator-pole-lap 
time of 191 ± 17 days. The result re ported here differs s lightly from 
the value of 199 ± 13 days given in Barnes et al ] <2005l) . but agrees 
within the uncertainties. The only difference here is that we have re- 
determined the optimal inclination angle, and use 70° rather than 
55°. 



7 DISCUSSION & CONCLUSION 

The images of Speedy Mic reveal it to be one of the most active 
stars to have been imaged to date. This is largely due to the high 
S/N ratio of the data combined with high spectral resolution of ~5° 
at the equator. Unfortunately S/N variations of the data and a lower 
degree of phase coverage on July 23 make comparison of some 
features difficult or impossible. In Fig.|5]we have reconstructed an 
image of Speedy Mic using all nights of data and the optimised 
differential rotation shear. 

7.1 Spot distribution and evolution 

In order to obtain a better impression of the surface distribution 
of spots, we have also plotted the projected disc of Speedy Mic 
at four key phases. This gives a clearer representation of the polar 
regions which are severely distorted in the Mercator projections. 
The plots make it clear that there is in fact no single uniform po- 
lar spot, as i s often seen on late-type rapidly rotating stars (e.g. see 
assmeier 2002). Except for the large de-centred spot/spot-group 
at phase 0.0 - 0.25, the high latitudes and the polar regions appear 
to comprise a number of resolved closely packed spots. A similar 
de-cent red spot was found on the young a Persei cluster G dwarf, 
He 520 1 B arnes et alll998lE00l[) . Whether a stellar magnetic cycle 
is responsible for the change in the morphology of spots in the po- 
lar region remains highly speculative. If meridional flows however 
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Figure 5. Image reconstruction using 2002 July If 
0.25 phase intervals. 



19 & 23 combined. Differential rotation included. Top: Mercator projection, Bottom: Projected disc at 



are responsible for carrying erupted magnetic flux to high latitudes 
fMackav et alj2004h . a persistent active region at a given longitude 
may be responsible for the appearance of a decentred spot. 

The images are remarkably similar to the images presented in 
B01. Observations were much more sparse during the 1998 obser- 
vations, and were made using more than one telescope and instru- 
mental setup. Nevertheless, the uniform distribution of spots at all 
latitudes is found in both sets of observations. The increased filling 
at high latitudes in the 2002 data may be partly real, but may also be 
the consequence of using a slightly higher inclination for the image 
reconstructions (see B01). 

Despite its rapid rotation, examination of the images in Fig. [3] 
shows that many features are repeated on both one and five day 
timescales. The spots may erupt in groups which are centred at 
phases of approximately 0.1, 0.3 and 0.5. Between phases 0.5 and 
1.0, such grouping is less clear. We caution that one should not 
be tempted to think that the smeared arcs (see ij6j in Fig. |4] con- 
nect spots together in a group. These arcs are regions of greatest 
uncertainty, and probably a consequence of the maximum entropy 
reconstruction process. 

Since observations of Speedy Mic were made by WSW05 be- 
tween two and three weeks after those presented here, we are in a 
position to examine evolution of surface features on this timescale. 
The observations presented by WSW05 were made with the UV- 
Visual Echelle Spectrograph at the 8-m Very Large Telescope, 
Kueyen. Deconvolution over a single spectral order enabled sin- 
gle profiles with S/N of several hundred to be obtained, and the 
images were derived using the updated CLEAN-like algorithm, 
first applied to obs ervations of the K0 dwarf AB Dor iKursted 
1 19931 Ktirst er et alJll994l) . The deconvolution technique applied 
by WSW05 differs from that used in this publication in that a tem- 



plate spectrum is iteratively convolved with an adjustable rotation 
profile, or broadening function, until an acceptable fit to the ob- 
served spectrum is obtained. This potentially enables a better fit 
to the observed spectrum to be derived as compared with the least 
squares deconvolution technique. The template spectrum method 
may be more efficient in boosting the effective S/N ratio than the 
least squares method, but the deformed rotation profile will be ef- 
fectively the same in both instances. WSW05 applied this method 
to two independent single echelle orders, yielding profiles with S/N 
ratios of 300 - 500. 



Whereas maximum-entropy regularisation is a continuous op- 
timisation problem, giving the smoothest possible image (i.e. the 
one containing the least amount of information), the CLEAN-like 
approach attempts a discrete approximation of the line profile de- 
formation (i.e. the starspots). Hence discrete filling factors are re- 
constructed, and it is argued that penumbral regions are more ac- 
curately represented using this method. Kurstei 1 1993) showed that 
CLEAN-like and maximum entropy reconstructions produce simi- 
lar maps, but that the appearance is more 'ragged' using the former 
method. This certainly appears to be the case when we compare 
our images with those of WSW05. However, it may be difficult to 
compare the two sets of images on such a small spatial scale since 
the S/N ratio of our data after deconvolution is nearly an order of 
magnitude greater as a result of using a larger spectral range. At 
S/N ratios of several thousand, the data should constrain the im- 
ages very well with the image entropy playing a smaller role in the 
reconstruction. As already discussed, the smoothing in our images 
must be largely due to the maximum entropy regularisation. It thus 
seems probable that we are not making the most of the high qual- 
ity deconvolved profiles. A re-examination of regularisation tech- 
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niques where high S/N ratio data are involved may be appropriate, 
but is beyond the scope of this paper. 

We find a remarkable correlation between the spot groups 
found in the present images and the image derived for 2002 Au- 
gust 02 by WSW05 (see Fig. 6, 6400A image). We suggest that 
WSW05 find the same spot groups at longitudes of approximately 
80°, 150°, 220°, & 280°. Since our data are phased on a different 
epoch (see f|6j to those of WSW05, we find that our phase 0.0 cor- 
responds to phase 0.24 in the WSW05 images. The longitudes of 
the spot groups in WSW05 correspond to phases 0.55, 0.34, 0.15 
and 0.98 on our maps, and are remarkably close to the locations 
described above, of 0.5, 0.3 and 0.1. The group at 0.98 may corre- 
spond to the group found at phase ~ 0.9 in this work. 

By 2005 August 7 (WSW05), the spot groups have evolved 
beyond recognition, apart from the phase 0.34 group which has be- 
come more extensive. The images presented here and in WSW05 
allow us to place a lower limit on the lifetimes of individual spots at 
five days. By contrast, most of the larger spot groups persist for at 
least two weeks on Speedy Mic, with only one group being present 
over a three week period. Unfortunately we cannot make any esti- 
mate of the real duration of spot groups without a longer span of 
observations. 

The two week timebase of the the data sets enables us to 
make an approximate verification of the differential rotation de- 
rived in this paper. The shift in the longitude position of spot groups 
is ~0.05 (18°), so the same features are appearing slightly later 
in phase by August 02. If we make the approximation that the 
spot groups are centred at latitude 45°, then the differential rota- 
tion derived in £16.21 implies Af2 = 0.016, giving a phase shift of 
0.0026 day -1 . Over the 15 day span from July 18 - August 2, we 
therefore expect a total phase shift of 0.039, which is in agreement 
with our estimate of 0.05. 

7.2 Comparison with other objects 

Other stars with similar sp ectral type such as AB Dor and L Q Hya 
(K0V; iDonati et al J 120031) . PZ Tel (KlV:lBarnes et aljEoOOI) and 
LO Pee (K5V: Lister et alj[l999l: iBarnes et alJl2005h . reveal very 
similar distributions of spots. Generally these stars exhibit either 
several large spots (e.g. LQ Hya and LO Peg) or several spot groups 
distributed relatively uniformly in longitude and latitude. In addi- 
tion to S/N ratio variations, it is clear that we must take into account 
the projected rotation velocities for each object when considering 
the degree of surface detail recovered. For LQ Hya, v sin i = 26 
km s" 1 and for LO Peg, v sini = 65 km s _1 (at approximately 2/3 
the instrumental resolution of the other objects). Both these objects 
exhibit large spots which no doubt comprise smaller unresolved 
features. Progressively more detail is resolved as we increase 
the projected rotation velocity, from PZ Tel (v sin i= 68km s -1 ), 
through AB Dor (vs'mi = 90km s" 1 ), to Speedy Mic (usini = 
132km s" 1 ). 

This resolution issue will also affect how much structure is re- 
solved in the poles and whether the polar spot appears solid, since 
the spatial resolution at high latitudes is lower. It is difficult to see 
a trend however since the two most extreme cases in the sample 
(in terms of rotation rate) possess the smallest (LQ Hya) or least 
extensive (Speedy Mic) polar caps. Given the long term nature of 
observation s for AB Dor, and the absence of polar spot in 1988 
lKursteretalll994l) . it seems most likely that magnetic cyclic ac- 
tivity affects the size and presence of a polar spot. It is difficult 
however to determine whether Speedy Mic exhibits cyclic activity, 
and whether the high latitude structure has changed significantly 



since the 1998 observations. Reconstructing the image with an ax- 
ial inclination of 70° results in an increase in the strength of the 
high latitude features. An effective (fractional) spot area of ~ 0.05 
at low-intermediate latitudes and ~ 0.15 at high latitudes are found 
(Fig.|3| right panels). If the present images are reconstructed at 55° 
it is found that the effective spot area as a function of latitude varies 
little, peaking at around 0.075. The images then closely resemble 
those found in B01. 

7.3 Prospects 

It is difficult to state whether or not a shift in latitude of the starspot 
groups has occurred over the timebase of observations. Simula- 
tions which focus on the role of enhanced magnetic flux emer- 
gence on the subsequent transport of flux at the surface indicate that 
the m ixed polarities seen on the K0 dwarf, AB Dor I Donati et al. 

12003), can be explained by meridional flows l Schriiver & Title 
200 lilMackav et all2004l) . If the 100 m s" 1 flows implied by these 
simulations are to be tested, daily observations over several days 
must be carried out. By investigating the possible change in the rate 
at which spots pass through the extent of the rotation profile in data 
space, it may be possible to obtain an estimate of the meridional 
flow patterns. Enhanced meridional flows may also be responsible 
for the presence of the polar cap. 

Unfortunately Speedy Mic is too faint for spectropolarimetric 
observations which may allow us to obtain some indication of the 
nature of the polarities at high latitudes. Since the contrast is so high 
in the optical between spot and photosphere, we would not expect 
to obtain any information from the spots themselves. We must turn 
to the next generation of infrared spectrographs which will allow 
us to use lines which only form in the spots to measure the polar- 
ities of the spots themselves. In the near infrared at ~ 1.7/im, the 
blackbody flux ratio is of order 2, whereas it is closer to 12 in the 
optical. 
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